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Модель многостадийного заболевания 
для задач медицинского страхования. 


У роботі запропоновано модель захворювання як багатостадійний компартментний процес. Використовується 
підхід на основі оцінки часу перебування особи (пацієнта) на кожній зі стадій. Розглянуто основні 
розподіли, пов'язані з часом перебування пацієнта на певних стадія захворювання - експоненціальний, 
Вейбула, Гомперца, Гомперца - Мейкхама. 
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В работе предложена модель заболевания как многостадийньй компартментньгй процесс. Используєтся 
подход на основе оценки времени пребьвания лица (пациента) на каждой из стадий. Рассмотреньт 
основнью распределения, связанньте со временем пребьтвание пациента на определенньтх стадия заболевания -- 
зкспоненциальноге, Вейбулла, Гомперца, Гомперца-Мейкхама. 

Ключевьге слова: математические модели, страховая медицина, компартментньтй процесс. 


Вступ 


Розрахунок актуарних показників у задачах страхової медицини вимагає роз- 
робки конструктивних математичних моделей багатостадійних захворювань. Переважно 
модель розглядається як (т -- 3)-стадійний компартментний процес. Припускається, 
що особи знаходяться на одній із цих стадій. Усі вони можуть бути вразливими (ста- 
дія 0); вони можуть бути на одній з т стадій захворювання (стадії ), ) - І, т ); вони 
можуть бути невиліковно хворими (або ж стадія повного прогресування захворювання, 
як наприклад СНІД або стадія раку ТМ-І - стадія т--1), або ж померти (стадія 
т- 2). Для різних видів захворювань стадії можуть класифікуватись по-різному. Так 
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для інфекційних захворювань початкові стадії пов'язані з інфекційними та латентними 
періодами, для онкологічних - це ранні стадії раку і т.ін. Особи можуть переходити 


зі стадії ) нанаступну стадію )--1, ) -0, т -- І з певним перехідними ймовірностями. 
У будь-який момент особа на стадії ) може померти з причин, не пов'язаних із за- 
хворюванням, тобто вона може перемістися зі стадії / безпосередньо на стадію смерті 


з перехідною ймовірністю и,, ) -0, т --1. Схема моделі наведена на рис. І. 


У 
Стадія 0 Стадія 1 Стадія т 1 


Стадія т- 2 


Рисунок 1- Компартментна модель захворювання 


Існує два основних підходи до побудови математичної моделі на основі ком- 
партментної моделі на рис. 1. Перший підхід розглядає чисельності індивідуумів на 
кожній стадії /. При цьому вони можуть бути описані як апаратом диференціальних 
рівнянь, так і ланцюгів Маркова з одиничним приростом (процес народження-смерт- 
ності) (1-31. Розв'язок таких диференціальних рівнянь та побудова реалізації ланцюгів 
Маркова становлять певні обчислювальні труднощі і не дають жодних аналітичних 
результатів. 

Другий підхід пов'язаний з розглядом часу очікування особи на кожній стадії. 
Отже, позначимо через У, час, який особа перебуває на стадії і до тих пір, поки вона 


не перейде на стадію (1 --1),і «0, т--1. Поряд з цим нехай (, - час, який особа пере- 
буває на стадії і до настання смерті, тобто до моменту коли особа переходить зі 
стадії і безпосередньо на стадію смерті т --2. Припускають, що величини (/, і У, є 


незалежними. 
Введемо наступні випадкові величини: 


Н, « пип(,, У), із 0, т - фактичний час перебування на стадії і та величини 
У «П,-У,. 

Очевидно, що якщо У » 0, то особа переходить із стадії і на стадію і-І. 
Якщо ж У/, « 0, то особа на стадії і помирає, тобто переходить на стадію т-2. 

Позначимо через У, величину загального часу очікування для особи, яка була 


на стадії і готова, щоб перейти на стадію 7-1. Тобто нам потрібен загальний час, про- 
ведений на стадіях 1і,...., / до моменту залишення стадії ) , щоб перейти на стадію ()--1). 
Тоді маємо: 


узУНн, теіздтяі. 


) 
каі 


-. 
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Позначимо через Х(т) є 10,1,... т-- 2) - стадія, на якій індивідуум перебуває в 
момент г. Метою роботи є знаходження ймовірностей: 
4: (1) - РІХ() - // Х(0) -ірприт»0,де іх ), і,/ -0,т-2. 


Основна частина 


Теорема 1. Перехідні ймовірності 4, (1), і 2 ), і, - От -- 1 можуть бути розра- 


ховані за співвідношенням: 
45 (1) - Р(5(1) - 1/5 (0) - і) - РІЙ, » 0,К -і,)-Цх (РУ, 0 ка, -1)- 
РЕА РТ 0 А); 


Доведення. Подія ГУ; » ГУ, » 0, К «і, - 1) відбувається у двох випадках: 
1) особа знаходиться на стадії ) в момент часу 1, тобто 5(1) - ), або 
2) особа знаходиться на стадії ()-Ї) в момент ї, причому загальний час 


очікування У, : 


7-12 5, і особа обов'язково перейде на стадію ) дещо пізніше. 


Звідси: 


РУ, з, »0,каїі,) -1) а РО/5П)- ШОУ, р», 0, ка і,)-Ще 
- РІ5() з ПУ, »0,кеі, | «ЛРУ ЗМ, 0, кі, і. 
Використавши формулу умовної ймовірності, після перегрупування маємо звідси: 
Р(5() з) і У» 0, ке) Пе РМ, з 0, ке і, -1х (РУ, з МУ, 0, а, і -1)- 
-РІУ р, з 0, сі, -Щ). 


Далі звернемо увагу, що подія (5(1)- //5(0)- і) може бути переписана як подія 


(5 (17) - ) 1 особа «стартувала» зі стадії і, вижила і успішно дійшла до стадії 7). Тобто: 


(51) « //5(0) «де Р(5(1)-«/ і М, »0,Каїі,) - 1). 
Теорему доведено. 
У подальшому в роботі буде зроблено спробу отримання конструктивних виразів 
для розрахунку перехідних ймовірностей 4, на основі Теореми І. 


Спрощена чотиристадійна модель 


Розглянемо спрощену модель у випадку тп - 1, яку графічно представлено на рис. 2: 

Компартменти моделі мають такі біологічні трактування. Стадія () представляє 
здорових, але вразливих осіб даної популяції (в епідеміологічних моделях - компарт- 
мент 5 ). Стадія 1 - стан первинних стадій захворювання (в епідеміології - латентний 
інкубаційний період - компартмент І,, ВІЛ, в онкології - рання стадії раку). Стадія 2 - 
стан розвитку і прогресування захворювання (в епідеміології - симптоматичне або 
асимптоматичне інфікування, СНІД, в онкології - пізні стадії раку). 

Зауважимо, що модель на рис. 2 може бути наближена до реальних за рахунок вве- 
дення додаткових компартментів (наприклад для усіх стадій раку 71,72, ТМ3, ТМ 4 ). 
Надалі обмежимося застосуванням моделі для опису онкологічних захворювань. Для 
цього зробимо такі припущення щодо розподілу величин Й, 1 У: 
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нні поні | 


«Р. 


Рисунок 2 - Чотиристадійна модель захворювання 


Припущення 1. Вважаємо, що величини (/,, і- 0,2 мають розподіл Гомперца -- 
Мейкхама СМ (и, |, й,» П,), де параметр и,, представляє незалежну від віку смерт- 


ність від захворювання, у, ,, П; - параметри масштабу та форми в розподілі смерт- 


ності через вік. Тобто функції щільності розподілу Ї,, і - 0,2 мають вигляд: 


з Да 
Ло з (шо Ж Мі2е" уехр| - Ма Я п. - о : (1) 


Припущення 2. Вважаємо, що величини У,, У, відповідають закону розподілу 
Вейбула Уеї(а.,, В,), і - 0,1, де параметр а, - параметр шкали, В, - параметр форми. 
Відповідні функції щільності розподілу мають вигляд: 

у а ваРч а уехр І (п/о) і 2 04. (2) 

Такий вибір закону розподілу для У,, і - 0,1 грунтується на тому, що закон Вей- 
булла вже використовувався у моделях багатостадійних процесів захворювань, на- 
приклад (1-3). З метою отримання представлення для 4,(г) нам потрібно здійснити 
ряд допоміжних викладок. 

Розподіли Н,. Випадок і - 0,1. З означення Н, та незалежності Ї, та У, випли- 
ває співвідношення: 


РІН, з 1) є Р(ті (0, У у зтре РМ, з, У, зт) аю - шия ро г х 
(3) 
В; 1- В; із та 
ХР сла Ре ехрі| - мого (од НО -е г» 0. 
Звідси функція щільності розподілу для Н., і - 0,1 становить: 
р 4 це В; Вал пи Ре і и, 2 пи 
Лу ( езяоРІН, збе (пор) те" ехрі з мі - (1/0) --9019671Ї; (4) 
і а бод; ; 
1»0. 
Розподіли М/,, і - 0,1. У подальших викладках нам необхідно отримати значення: 
РІ, » 0) - Р |особа перейшла зі стадії і в стадію (і --1) /. 
Тобто: 
РІ, » 02 | учи, б, (5) 
0 
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де Ї у (1) - функція щільності розподілу величини ГП, - У, яка може бути от- 


римана на основі теореми про згортку |41: 
ФЕЇ, баз 0-97, Оз). (6) 
Звідси, підставляючи в (6) вирази (1) 1 (2), маємо: 


рн повно аох 


РО Гама ЧИ2е") "В ра Я (7) 


і 


хехр- ((з - 1// а)? |аз 
Звідси: 


РІ, » 0) - Гоа Ж ою - о 
01 


празету(рівоюминЙх ||) 


хехр- (з -гу/ а) дз, з є 04. 
Зазначимо, що в загальному випадку (8) - інтеграл, що не береться. Для конкретних 
задач він може бути обчислений чисельними методами. Розглянемо приклад, коли: 


шаг20,В,-ї, (9) 
РР (10) 
Ки, 


де К - 1,2, 3... є невід'ємне ціле. 
Умови (9) обмежують вашу модель до випадку, коли У, експоненціально роз- 


і 
поділена з параметром а, а У, розподілена згідно з розподілом Гомперца з параме- 
трами и,» та 7,. Умова (10) вимагає, щоб параметр шкали о. вимірювався в одиницях 
параметра форм розподілу Гомперца т, 1 навпаки. У такому випадку для обчислення 
щільності розподілу У/, може бути отримана конструктивна формула. Наприклад 


коли К - 3, то: 


2. т; 
12 зе 
( т; - о т; " Е; И, М 2 ) 
Ні, 2 се Р ) і 
и РАК т; 2 РАНІ у 
и3 ее ауд, ее 
(0) - ша З рах х 
5-200 ЦІР 
Ні, 2 
реа 3 е'"їїу, зб'яЖтИі 
хире'елтич Го п? ролів" ж ЕЙ, Мі 2 унде Міо Пі ух (11) 
Ні 2 тріо Мі" 
е РА БА Р 
ПТ; 

У другому варіанті чотиристадійної моделі замість Припущення 2 припустимо 
наступне. 

Припущення 2. Вважаємо, що величини У,, У, відповідають закону і на 
Гомперца - Мейкхама СМ (с. |, Оз, В, ). Тобто функції щільності розподілу У,, і - 0,1 
мають вигляд: 

а. 
о о а за пр 12) 


164 «Искусственньєй интеллект» 172012 


Модель багатостадійного захворювання для задач медичного страхування ам 


Розподіли Н,. 
РІН, »1ре РИ ті (У, У ) 2 1) - РІО, 21,У, ачек- най оту Роз а- го 
т; 
1»0 


і і і 


о, о, 
х "ЕХ нАІ- 2 че ГОРА о КА Ва зелунанІ- що» 
Р Р 
Звідси функція щільності розподілу для Н,, і - 0,1 становитиме: 


ЛФ) ЗУ РІН, зафе(ЧИ, зол тує зВе""ух 


(13) 
Йо лях, Обі2 Ви 
мехр- і НО Рено НВО «РЕНО 
т; В; 
Розподіл У/, і «0,1. 
Введемо допоміжні РипаЮнеІ величини: 
шов ззаца казйозуг, 14) 
При цьому: 
РІЙ з 0) - Р(е"" зе" ре Ру, » 0). (15) 
Застосовуючи формулу для обчислення щільності розподілу випадкової вели- 


чини, яка є монотонною функцією від іншої випадкової величини (4), маємо: 


Ра а-т) ря 


чт Міо Мі 


ері оте пі - 
ее з(фажиої Й е Є 


1 1 
Уа0- гм Пп) з вна г Шо" у - Па 


ті "Ні ціа МП; 
31 ,1 е 2 


1 
Б рано ЖИ 


р а , зо іосуфії 
ел (1) «б ва і ОР Ї ча Ві І | 7) 


Застосовуючи теорему про згортку, маємо: 


пе з) у аа Ду (о 8), (здав з 


т і "Но і2! фе -зуі 
-Ї ни т наб з) Хі 8 ее з: (18) 


хацунард УС 1 ебнаЙотСЯ дб, 
Розглянемо частковий випадок (12) з параметрами: 
ша 70,7 0,п, еВ, 22, 19) 
що відповідає РР Гомперца з параметром форми, рівним 2. У такому випадку: 


1 
НО ай НО 2) 1 І за у з) 
Лу, зу п) пев 0 М ЛІ ех "їв Мао о МД е? з 


знака) (20) 
А Шо 2 ен - з 2452р-2л ен( -- дує и 
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У третьому варіанті чотирастадійної моделі замість Припущення | 12 припу- 
стимо таке. 


Припущення 1. Вважаємо, що величини 0, - 0,2 мають експоненціальний 
розподіл Ехр(и,), тобто функції щільності мають вигляд: 
АФеце"", 1»0, із-02. (21) 
Припущення 2" збігається з припущенням 2! (тобто У,, У, розподілені відповідно 
до закону Гомперца - Мейкхама СМ (и, й, П;), 12 0,1.) 


Розподіл Н.. 


зії ОС; ї 
РІНО зр РІШИОГ У уз Пер ЗУ зароє "БУ б -ей Б 


2 Ва 
чені ана оза і ) 1»0. 


і 


Звідси функція щільності розподілу для Н., і - 0,1 становитиме: 


д за 
тамооує зовн »де(ащ о це? )х 
(22) 


ке - (си ..; а-еїу кає 0. 


і і 


Розподіл Й. Використаємо допоміжні випадкові величини (14). Застосовуючи (15), 


маємо: 
АД (з у (піг) з зр нен зн; 
о , го шВі 
РРО С зд оо а (23) 
Застосовуючи теорему про згортку, маємо: 
-ші-1 
лате | 7.55) з | шоу х 
(24) 
х (ау Оз Снзу Сн пе" РетОо до, 
Покладемо а, - 0 (тобто розглянемо лише розподіл Гомперца). Маємо: 
а зії Вічної 
-5йі -ші-! Ві В; 
| ета, и; зе і 
2 (1) б пп зу ні - 
5-00 В ; В й 
М (25) 
ИОС, 2 з 
ш- 8 . Ж 2 
Звідси для як завгодно малого є » 0: 
аа аа 
В; оо -іРі Чі В; -є Іі 
Хе е (-1) і а,е"" е 
і є Ї і Ї 
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Модель багатостадійного захворювання для задач медичного страхування ам 


З метою спрощення технічних викладок введемо спеціальну функцію. 
Означення 1. Інтегральною показниковою функцією порядків а та Р називається: 


в аб) з | у 07) 
-о 5 
При г20 інтеграл в (27) розуміється в сенсі головного значення, тобто: 


ва ій пре 
е е 
Еї(1,а,Б) - йт Ї --у4У Я йт ГУ з (28) 
єзно 7, 5 есзжОГ 5 
Зауважимо, що існує зв'язок Е(1,а, Б) з інтегральною показниковою функцією 
Ейлера а саме: 


ЕІ) « ЕМо,1,1). 


Отже, використавши Означення | з (26), отримуємо: 


О 2 


з НОС; зе 


В; 
РІ зер е (Тв ЕСЕ Вон 1. (29) 


і 


Висновки 


Отже, в роботі запропоновано модель захворювання як багатостадійний ком- 
партментний процес. При цьому, на відміну від традиційного підходу в медичному 
страхуванні, який розглядає ймовірну чисельність осіб на кожній стадії і задовольняє 
задачі епідеміології, тут використовується підхід на основі оцінки часу перебування 
особи (пацієнта) на кожній зі стадій. Такий підхід є загальнішим, оскільки дозволяє 
охопити ширший клас соматичних захворювань. Також у роботі розглянуто основні 
розподіли, пов'язані з часом перебування пацієнта на певних стадія захворювання -- 
експоненціальний, Вейбула, Гомперца, Гомперца -- Мейкхама. Слід зазначити, що роз- 
рахунок перехідних ймовірностей процесу багатостадійного захворювання в класі таких 
розподілів часу перебування можливий лише з використанням спеціальних функцій. 
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ам Марценюк В.П., Климук Н.Я. 


У.Р. Маті5епуик, М.Уа. Кіутик 
Моає! о) Митуаєе Різеазе ог те Ргобіету ої Меаїсаї! Ппзигапсе 


Тре таїп їдеа ої Фіяє 5(иду 15 а плоді! ої 4їзеа5е а5 а плиїшякаєє соптрагітепі 
ргосез58. ТРре иц5ед арргоасі 15 Ба5ед оп йте аз5ез5тепі 5репі Бу а рег5оп (райепі) аг еаср 
зкаєе. Тре пзаїп Фі5пібибоп5 аз5осіатед улір Ше те Ше райепі 5їауз аї сегіаїп 5(аєе5 ої 
ре дїзеазе, 1.е. ехропепиаї!, Метиі, Сотрегі?, Сотрегі7-Макерата, аге 5їидтеа. 

Треге аге їмо Базіс арргоаспе8 іо Бийадтпє плайептайса! піодеїя Ба5ей оп соппрагі- 
тепі плодеі. ТБе бг5і арргоасП соп5ідег5 Ше питкіег ої іпдутдцаїз аї гасі 5(аєе. Три8, Шеу 
сап Бе дезсгібед а8 аррагашя ої іНегепиа! едиайопя апа МагКоу сраїпя упр а ппік іпсгеазе 
(річШ-дсай ргосез5). ТПе 5о0Іипоп ої зисп ФїНегепіа! едиайопя апа Ше соп5ігисйоп ої ре 
ітріетепіайоп ої Магкоу сраїп5 плаке сегіаїп согарикайопа! дїйбсиїйеє ап 9іуе по 
апатубса! гезиїйя. Те зесопа арргоасі 15 сопсегпеай ул ре сопзідегацоп ої маїйпо, йте8 ої 
а рег5оп аг еасп 5каєе. Тре аїпа ої ія 5(иду 15 (о ргоробе аїсогійптя їог (ре саїсціайоп ої 
іап5ійоп ргобабіййез, пзіпє Ше арргоасіп Ба5ед оп ууаїйпя те. ТРі5 арргоасі 15 плоге 
сопатоп, Бесацзе її аПом/з соуегіпе а у/1дег сіав8 ої 50тайс ді5еа5ез. 


Стаття надійшла до редакції 24.11.2011 І. 
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